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^ ; ABSTRACT 

Aims. To investigate structure of the hot object in the symbiotic prototype Z And during its major 2000-03 active phase. 

Methods. Analysis of the far ultraviolet, optical low- and high-resolution spectroscopy and UBVR photometry. Reconstruction of the spectral 

energy distribution (SED) during the outburst. The Raman scattering process. 

Results. At the initial stages of the outburst the hot object was characterized by the two-temperature spectrum (a warm stellar radiation and a 
strong nebular emission) with signatures of a mass-outflow at moderate (~ 100-r200kms"') and very high (a; 1000-^2000 km s"') velocities. 
The corresponding structure of the hot object consists of an optically thick, slowly-expanding disk-like material encompassing the accretor at 



Oh 
I 



O the orbital plane and a fast optically thin wind over the remainder of the star. The disk-like shell persisted around the central star until 2002 



August as was indicated by the eclipse effect. Then a significant dilution of the optically thick material and evolution of a fast wind from the 



CZ2 hot star, concentrated more at the orbital plane, were detected. A striking similarity of [Fe vii] A60S7 and Raman A6S25 profiles at/after the 
dilution of the disk suggested their origin within the interaction zone where the winds from the binary components collide. 
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1. Introduction mated a roughly constant luminosity and a factor of 10 lower 

temperature, T^, for the hot object during eruptions with re- 

ZAnd is considered as a prototype of the class of symbiotic . . .... ^ ■ » . inoc 

, ,. r , .,.^TTT • . spect to quantities from quiescence . As co ncerns to the 1985 

stars. The binary composes of a late-type, M4.5 III, giant and a . trr: — rr ; — r—r 1 i,„„ j . j j 

■' . ^. , , . , , -T^o , eruption. IMikolaiewska & KenvonI ( 1199a) supported the de- 

white d warf accreting from the giant s wind on the 758-day or- • ^ i u j i- en j ^■ 

, . , rr — ; "„ — rrrrnTirrr; — ' „ — ' 1 crease in Th also by a decune of Raman scattered Ovi hnes 

bit (e.g. lNussbaumer & Vogel ll989tlMikolaiewska & Kenvon I . ^^otc . nnoo t-u- .u » » j- j • j » -i u 

I " V l »T — r TTZ — 77^ at A6825 and A/088. This out burst was studied in detail by 

119961) . More than 100 years of monitoring Z And mrst records r: — — ~ — m Jinnee r> c i. • i . 

— ^ ■' " IFernandez-Castro et al. I r 1995). By means of the ultraviolet 

from 1887) demonstrated an eruptive character of Its light curve . j j - u ■ j - . j 

' r b spectroscopy and radio observations, they indicated a con- 

(LC). It displays several active phases, during which fluctua- . , , , , . .... , ^ 

^ ±- ' o siderably cooler UV continuum with an additional attenua- 
tions ranges in amplitude from a few tenths of a magnitude . , , nrvrv a • j i. j v • »u j- 
. . jr-i u u /inic '■1°'^ below 1300 A accompanied by a decline in the radio 



to about 3 magnitudes. During the strongest outbursts (1915 
1940, 1960, 1985, 2000) the star's brightness reached values 



flux during the 1985 optical maximum. In the line spectrum 
they measured a broadening of all emission lines with signa- 



around m^„ ~9 and Wvis ~ 8.5 mag with the recurrence time ^ c r, ■ ^ r^ir^ i- 

^'^ tures of P-Cygni type o f pronles for strong resonance lines. 

of 15 to 25 years (see Fig. 1). Maior outbursts are charac- 17 — \ ' , I ,innc> • . . j .u f . 

, . , ' . , , IFernandez-Castro et al. I (11225) interpreted these features as a 

terized by a decrease in the eltective temperature of the hot i i • r n i ■ i i n . j i 

, , r. r • • result of the ejection of an optically thick shell at moderately 

obiect and the fluxes of some high ionization emission lines. , . , , r^rvrv orvrvi -i r. j- i. i 

f 1 high velocity of 200-300 km s . Radio observations also pro- 
Analyzing previous large eruptions, pvIikolaiewska&Kenvom • . . • . t t . 

J ' i ^ o ±- ' f 1 — I vided evidences of a mass-outflow from the active star At the 

(119961) found that an A-F absorption spectrum character- • i ^ f u »i, inoc j onnn . a- a 

' ^ ^ ^ initial stage of both 1985 and 2000 outbursts the radio flux 

ized the 1940 optical maximum and during the 1960 out- . j u i if ■ a ^ 

^ ° dropped below values from quiescence, then recovered along 
burst the system resembled an FO-2 supergiant. They esti- ^, i j v a u a ^ . .u • .11 
^ ■' the optical d ecline and reached close to the quiescent level 

* Visiting Astronomer: Capodimonte Astrophysical Observatory or above it jFernandez-Castro et all 1 19951 iBrocksoPP et al. I 
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|2004 . For the 1984-85 activitv. lFernandez-Castro et al. I ( ll995l) 

explained the anticorrelation between the visual and far-UV 
LCs by the extinction of UV radiation from the central source 
by a false photosphere produced by the ejection of a shell of 
dense material. They compared the subsequent rise of the radio 
emission with a model, in which the ionizing source is turned 
on instantaneously at the center and an ionization front propa- 
gates out thr ough the wind from the hot star For the 2000-03 
active phase. iBrocksopp et al. I(l2004l) revealed a jet-like exten- 
sion at 5 GHz on the 2001 September MERLIN image. They 
associated this transient feature to a blob of matter ejected at 
the beginning of the optical outburst, emitting most probably 
through thermal bremsstrahlung of ionized hydrogen gas. 

On 2 000 September 1st Z And entered the recent major 
outburst dSkopalet al.ll200d ). It reached a maximum bright- 
ness around the mid of 2000 December (U ~ 8.4, B ~ 9.3 and 
V ~ 8.8 mag) following a gradual decrease t o quiescent values 
at/after the mid of 2003 fek opal et al. |2004 >). In the summer of 
2002, Skopal (2003a) observed a narrow minimum in the LC 
at the position of the inferior conjunction of the giant. He asso- 
ciated this effect to the eclipse of the active object by the giant 
and determ ined the orbital inclination of the system to ^ 76°. 
ISokoloski et al. I ( I2005ID observed the 2000-03 outburst in the 
X-rays, far ultraviolet, optical and radio. They focused on the 
nature of the outburst, and described observational evidence for 
a new kind of the outburst that involves both an accretion-disk 
instability and an thermonuclear shell burning on the surface of 
an accreting white dwarf. 

In this paper we present and analyze our original 
spectroscopic and photometric observations complemented 
by the archival spectra obtained with the Far Ultraviolet 
Spectroscopic Explorer (FUS E) carried out during the 2000- 
03 major outburst of Z And. In Sect. 3 we analyze the data. We 
model the SED (Sect. 3.3) and examine the Raman scattering 
process (Sect. 3.4). In Sect. 4 we summarize main results to re- 
construct a basic structure of the hot object and the surrounding 
environment during the outburst. 

2. Observations and reductions 

2.1. Optical spectroscopy 

Our spectroscopic observations were taken during the recent 
major 2000-03 active phase at different observatories. 

Prior to the maximum of the star's brightness a 
low-dispersion spectroscopy was secured at the Special 
Astrophysical Observatory of the Russian Academy of 
Sciences (SAO in Table 1) with the long slit UAGS spectro- 
graph equipped with a fast Schmidt camera and CCD detector 
(530x580pixels of 18x24//m size) mounted at the Cassegrain 
focus of the Zeiss 1 -m telescope. Linear dispersion was 1 .5 and 
3.1 A pixel ' in the blue and red part of the spectrum, respec- 
tively. Basic reduction including extraction of ID-spectra from 
CCD images and wavelength cali bration were made by using 
the software-package of Vlasvuk] ( Il993l) . 

At the maximum of the star's brightness a high-dispersion 
spectroscopy was secured at the Asiago Astrophysical 
Observatory (OAC) with the REOSC Echelle Spectrograph 



Table 1. Log of spectroscopic observations 



Date 


Julian date 


Phase* 


Rej 


gion 


Observatory 




JD24... 




[nm] 




02/11/00 


51851.27 


0.143 


409.2 


- 494.8 


SAO 


04/11/00 


51853.34 


0.146 


413.2 


- 495.0 


SAO 


05/11/00 


51854.37 


0.147 


413.8 


- 495.0 


SAO 


05/11/00 


51854.39 


0.147 


482.0 


-563.1 


SAO 


16/11/00 


51865.00 


0.161 


90.5 ■ 


■ 118.8 


FUSE 


27/11/00 


51876.00 


0.176 


90.5 ■ 


■ 118.8 


FUSE 


10/12/00 


51889.41 


0.194 


459.5 


- 909.2 


Asiago 


11/12/00 


51890.35 


0.195 


459.5 


- 909.2 


Asiago 


12/12/00 


51891.30 


0.196 


459.5 


- 909.2 


Asiago 


15/12/00 


51894.00 


0.200 


90.5 ■ 


■ 118.8 


FUSE 


05/01/01 


51914.92 


0.227 


580.0 


- 700.0 


OAO 


08/01/01 


51917.90 


0.231 


450.0 


- 570.0 


OAO 


27/05/01 


52057.26 


0.415 


600.0 


- 700.0 


OAO 


22/07/01 


52113.30 


0.489 


90.5 - 


■ 118.8 


FUSE 


30/09/01 


52183.30 


0.582 


90.5 - 


■ 118.8 


FUSE 


10/02/02 


52315.93 


0.757 


580.0 


- 700.0 


OAO 


05/07/02 


52461.50 


0.949 


90.5 - 


■ 118.8 


FUSE 


07/08/02 


52494.25 


0.992 


580.0 


- 700.0 


OAO 


22/10/02 


52570.50 


0.093 


90.5 - 


■ 118.8 


FUSE 


02/02/03 


52673.45 


0.229 


580.0 


- 700.0 


OAO 


31/07/03 


52852.50 


0.465 


580.0 


- 700.0 


OAO 


04/08/03 


52856.80 


0.471 


90.5 - 


■ 118.8 


FUSE 


02/09/03 


52885.50 


0.509 


580.0 


- 700.0 


OAO 


13/11/03 


52957.43 


0.604 


625.8 


- 677.0 


Ondi^ejov 


18/12/03 


52992.28 


0.650 


625.8 


- 677.0 


Ondi^ejov 



y/)Min = 2414 625.2 + 757.5 x£ iskoDalll998h 



(RES) equipped with a Thompson THX31156 UV-coated CCD 
detector (1024x1024 pixels of 19 pm size) mounted at the 
Cassegrain focus of the 1.82-m telescope at Mt. Ekar. Two 
very different exposures of 120 and 3 600 s were applied to 
obtain well defined Ha profile and the underlying contin- 
uum. The RES orders were straightened through the software 
developed at the Astronomical Observatory of Capodimonte 
in Napoli. Thereafter, the spectroscopic data were processed 
by the MIDAS software package. Dispersion of the RES was 
1 A mm" ' (0. 1 9 A pixeL ' ) at Ha. The spectra obtained during 
each night cover the wavelength range of 459.5 - 909.2 nm. 

After the maximum and the following transition to qui- 
escence a high-dispersion spectroscopy was secured at the 
Okayama Astrophysical Observatory (O AO) with the H igh 
Dispersion Echelle Spectrograph (HIDES, fizumiurall 19991) at 
the f/29 coude focus of the 1.88-m telescope. Dimension of 
the CCD (EEV 42-80) was 2048x4096 pixels of 13.5//m size. 
The spectral resolving power was 50 000 or 68 000. The dis- 
persion was ~0.02ApixeL' (~ 1.2kms"')and ~0.03ApixeL^ 
(~ 1.4kms-') in the region of 450 - 570 and 580 - 700 nm, re- 
spectively. Bias and instrumental flats were also measured and 
taken into account in the usual way. The reduction and anal- 
ysis was performed with the IRAF-package software and our 
own code to reduce spectra obtained with HIDES. For example, 
the scattered light from the instrument was also removed. The 
Th-Ar comparison spectrum was fitted with a 3rd-order poly- 
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nomial function, which reproduced the wavelength calibration 
to better than 10"^ A. 

At the post-outburst stage (the end of 2003) a high- 
resolution spectroscopy between 6258 and 6770 A was taken 
at the Ondfejov Observatory using the coude spectrograph 
of 2-m reflector and the BROR CCD camera with the SlTe 
800x2030 pixels chip. The resolution power at the Ha region 
was 13 000. Basic treatment of the spectra was done by using 
the IRAF-package software. 

A correction for heliocentric velocity was applied to all 
spectra. The journal of spectroscopic observations is in Table 1 . 

2.1 .1 . Flux calibration of the spectra 

We converted arbitrary flux units of our spectra to fluxes in 
ergcm"^ s"' A"' with the aid of the simultaneous UBVR pho- 
tometry. Following points are relevant; (i) We summarized pho- 
tometric observations carried out (nearly-)simultaneously with 
the spectrum under consideration, (ii) We calculated correc- 
tions Am, which are due to emission lines in the spec trum, to 
get m agnitudes corresponding to the continuum (see ISkopal I 
L2003b, for details), (iii) Then we dereddened the continuum- 
magnit udes with £r- v = 0.30 by using the extinction curve of 
ICardelli et al. ' ( 1989) and converted to flu xes according to the 
calibration of Henden & Kaitchuckl ( ll982l) . (iv) Finally, we de- 
termined the star's continuum by fitting the UBVR fluxes with 
a 3-order polynomial function and scaled the measured spec- 
trum to this fit. 



2.2. Far ultraviolet spectroscopy 

The F US E instrument was desc ribed in detail bv lMoos et al. I 
( I2OOOI) and lSahnowetal.l(l2000t). A brief useful description of 
the spectrograph is given in lYoung et al. I ( l2005l) . 

The data were processed by the calibration pipeline version 
3.0.7, 3.0.8 and 2.4.1. We used the calibrated time-tag obser- 
vations (TTAG photon collecting mode) taken through the large 
aperture (LWRS). Before adding the flux from all exposures we 
applied an appropriate wavelength shift relative to one so to get 
the best overlapping of the absorption features. Then we co- 
added spectra of individual exposures and weighted them ac- 
cording to their exposure times. Finally, we binned the resulting 
spectrum within 0.025 A. In order to determine the wavelength 
scale of the spectra, first a correction for heliocentric velocity 
including that of the satellite was applied. Second we selected 
interstellar absorption lines for which wavelengths are known 
and fitted the corresponding observed features. For example, 
in the region of the LiFlA spectr a we used N 1/11003.372, 
Si 11^1020.699, A I ^1048.218 (e.g. 'Rogerson & Ewefl "l984' 
Fig. 3 here). An accuracy of such calibration is of + 0.05 A. The 
data retrieved from the FUSE archive were processed mostly 
by own codes (AS). 

Dereddening t he FUSE s pectra we found the extinction 
curve of Cardelli to be valid only in the interval 

from the long-wavelength end of the spectra to about 1 100 A 
(open circles in Figs. 6 and 7; /iB v - 0.30, /?v = 3.1). Beyond 
these wavelengths to 1 000 A dereddened fluxes enlarges by a 



factor of 1.5-^2 with respect to quantities required by our mod- 
eling the SED (crosses in Fig. 7). Note that the extinction curve 
in the far-ultraviolet (1 250-^1 000 A) as suggested by Cardelli 
et al. (1986, Fig. 5 there) represents only an extrapolation and 
thus is more uncertain. As a result we used their mean extinc- 
tion law only for /I > 1 100 A. The short-wavelength rest of the 
spectrum was rectified to our model SED (Fig. 6). 

2.3. Broad-band photometry 

Photometric observations were performed in the standard 
Johnson UBV system with the Rc filter of the Cousins system. 
We used a single-channel photoelectric photometers mounted 
in the Cassegrain foci of 0.6-m reflectors at the Skalnate Pleso 
and Stara Lesna observatories. We converted our Rq magnitude 
to that of Johnson system by using transformations of Bessellj 
(1983). ZAnd was measured with respect to the comparison 
SAO53150 (BDh-474192; V = 8.99, B-V ^ 0.41, U - B ^ 
0.14, V - R = 0.16). Usually, 1-hour cycle contained about 10 
to 20 individual differences between the target and the compar- 
ison. This approach reduced the internal uncertainty of such 
the means to ~0.01 mag in the V, B and R bands, and up to 
0.025 mag in the U band. Individual measurements in a f orm 
of the table were already published by Skop al et al. I (l2004l) . 

In this paper we dereddened our observations with E^-v - 
0.30 and scale d relevant parameters for a distance of 1 .5 kpc 
dSkopal I2OO5I and references therein). 

3. Analysis 

3.1. Photometric evolution 

Bottom panels of Fig. 1 show the UBVR LCs of ZAnd cov- 
ering its recent outburst. The rise to maximum was character- 
ized by three rapid increases in the star's brightness and two 
plateaus, best seen in the U band. The initial rise between 
2000 September 1st and 7th (Af/=1.9, AB=1.5, AV^l.l) was 
probably caused by an increase of the hot object luminosity. 
The U - B index became more blue with respect to observa- 
tions made prior to the outburst. Consequently, the profile of 
the optical continuum became more steep. Its nature - nebu- 
lar or stellar (or a combination of both) - is important to un- 
derstand the trigger mechanism of the outburst. The nebular 
nature of the optical light at very beginning of the outburst 
would require an increase in Th to get the necessary amount 
of ionizing photons, while the stellar nature of the optical radi- 
ation would be associated with a creation of a shell surround- 
ing the hot star implying thus a decrease in T^- Unfortunately, 
we do not have relevant spectroscopic observations from the 
initial stage of the activity to distinguish between these possi- 
bilities. Nonetheless a comparison of colour indices with t hose 
of a b lack body (they are well above, see Fig. 1 of Sko pal I 
l2003b ) signals a contribution from other sources. The follow- 
ing short-term decline in the U - B index at constant B - V 
and V - R indices suggested a decrease of the hot component 
temperature. Then during the 'totality' of the U - B minimum 
(~ yD 2 45 1 800 - 2 45 1 865) all the colour indices were practi- 
cally constant indicating thus a constant temperature. Therefore 
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Fig. 1. T op: The historical 1 895-2004 ph otographic/B-ba nd/visu al LC of ZAnd. It is com piled from the data publi s hed by 
Pavne-Ga poschkin ( 1946), i Romanol (U^O), Mjalkovsk^ lll977l) . Belyakina ( 1985, \l99A |Mikolajewska & Ke nvonI ll 19961) . 



Skopal et al. > ( 2002b. 2004) and smoothed visual AFOEV data. Middle: A part of the multicolour LC covering the recent two 
major outbursts (1984-86 and 2000-03). Bottom: The UBVR photometry from the major 2000-03 outburst. Left panel shows the 
overall evolution in the LC. The eclipse is denoted by E. The right panel shows the LC at the initial stage of the outburst - note 
its stepped profile on the rise to the maximum. Optical/FUV spectroscopic observations are denoted by vertical full/dashed bars. 
The arrows mark positions at which we reconstructed the observed SED. 



the second increase in the star's brightness, which occurred at 
constant indices (~ JD 2 45 1 820, Fig. 1), had to result from an 
increase in the luminosity at a relevant expansion of the shell. 
Finally, the third brightening from the mid-November to the 
maximum (i.e. the ascending branch of the U-B minimum) 
indicated an increase of the hot object temperature and its lu- 



minosity. Quantities of the corresponding colour indices were 
probably due to a composition of contributions from a nebula 
and stellar source. Our SED model from the maximum sup- 
ports this interpretation (Sect. 3.3.2, Fig. 7). It is of interest to 
note that similar behaviour in colour indices was also observed 
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Fig. 2. Evolution in selected lines, mainly from the optical maximum. Fluxes are in units of 10"'^ erg cm"^ s"' A"' and radial 
velocities in kms Small numbers in panels with more profiles represent their shift with respect to the level of the local 
continuum. The arrows denote absorption components. The right panel shows evolution in the He 1/15876 line (y-tics are by 100 
units). The continuum from the giant, as given by the SED (Fig. 7), is compared to clarify details in lines of Hj6, [O m] /15007 and 
He I ^6678. 



27/11/00 




1020 1025 1030 1035 1040 



Wavelength [A] 

Fig. 3. Evolution in the Ovi/11032 line (LiFlA spectra) in 
an appropriate logarithmic scale. Note that the P-Cygni pro- 
file with a strong absorption developed throughout the whole 
active phase. Positions of interstellar lines are denoted by i. 



at early stag es of classical novae, p rior to the maximum (e.g. 
Fig. 3 of Chochol & Pribulla Il997l) . 



From 2001 ZAnd entered the decline phase during which 
both colour indices and the star's brightness were evolving to- 
ward quantities of a quiescent phase at the end of 2003 (Fig. 1). 
In the summer of 2002 the gradual decline was inteiTupted by 
the eclipse effect. This event revealed a high orbital inclination 
and the profile of the minimum suggested a disk-like structure 
for the eclipsed object (see Sk opal 2003a. in detail). In addi- 
tion, a transient ~ 0.6 mag rebrightening lasting from the au- 
tumn 2002 to the spring 2003 was indicated in the LC (see 
Sect. 3.6 in more detail). 

Our reconstructed SEDs, we present in the section 3.3, 
should aid to understand better the photometric behaviour we 
discussed qualitatively above. 

3.2. Spectroscopic evolution 

Figures 2-^6 show main characteristics of our spectroscopic 
observations along the outburst, from its maximum in 2000 
December to its end in 2003. We comment on them as follows: 
(i) The P-Cygni profiles of Hi lines of the Paschen se- 
ries and He i lines were present mainly at initial stages of the 
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Fig. 4. Evolution in the Ha, Raman /(6825 and [Fe vii] 6087 line profiles along the outburst, from the 2000-maximum to 2003- 
quiescence. Left panel shows profiles of the Ha line core (y-tics are by 500 units), while the middle panels show the extended 
wings of the profiles in detail. On the right panels y-tics are by 20 units and the horizontal dashed fines represent the approximated 
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Fig. 5. Observed LiFlA spectra from 30/09/01 and 04/08/03. 
The level of the continuum significantly decreased while fluxes 
of the O VI doublet increased by a factor of ~ 10 on 04/08/03. 
Note also the appearance of the Ne vi lines. 



outburst, prior to and at the maximum. The main absorption 
component was located at about -90 km s"'. However, com- 
plex absorption features to 300 km s ' cou ld be recognized 

for i8 502, ^8545 and ^8665. In the far-UV ISokoloski et al. I 
( l2005b detected P-Cygni profiles in C iii and P v lines (Fig. 6 
here) with terminal velocities of « -200km s"' on the first 
FUSE observations. They also demonstrated evolution from 
absorption to emission profile of these lines as the outburst 
progressed to its maximum. On the other hand, the resonance 
line Ovi/11032 shows a strong absorption in the P-Cygni line 
profile throughout the whole active phase (Fig. 3). This sug- 
gests its origin in the stellar wind of the hot star During quies- 
cence the absorptioii component was not present in the profile 
dSchmid et al. Ill99^ . In Sect. 3.4 we examine this line in de- 
tail. 

(ii) The H/3 wings extended to about 1000km s"' or more. 
The main emission core showed an asymmetrical profile with 
the peak placed at about +30kms"' (the systemic velocity 
is -1.8 km s"'). The asymmetry can result from superposition 
of a violet shifted absorption component created in the neu- 
tral gas at the front of the expanding material. The Ha pro- 
file exhibited similar behaviour. However, its wings extended 
to more than +2 000km s"^ and the asymmetry was more pro- 
nounced. At/after the spectroscopic conjunction the absorption 
component created a cut-out on the violet side of the wing 
at about -40km s"' (07/08/02) and -50km s"' (02/02/03). 
During the low-amplitude brightening (2002/03) we measured 
an additional extension of the Ha wings on our spectrum from 
02/02/03. We ascribe their origin to the stellar wind from the 
hot star (Sect. 3.5). 

(iii) At the maximum, the main emission core of the 
He II /14686 was shifted by about -75 km s"' at the peak and by 
about -105 km s"' at its bottom with respect to the reference 
wavelength. The [Fe vii] /J6087 line displayed similar charac- 
teristics. Its profile was very low and broad with the main emis- 



sion located around -80kms ' (Fig. 4). On the transition from 
the optical maximum to quiescence these lines became stronger 
being located at/around their reference wavelengths (see Fig. 4 
for [Fe vii] and Tomov et al. 2005 for the Hen line). The evo- 
lution in these lines is discussed in Sect. 4. 

(iv) He: /15876, /16678 and A1065 emission lines can b e 
used to diagnose the symbiotic nebula (Proga et al.1 Il994i) . 
We observed all these lines just at the maximum, but later, 
along the investigated active phase, only /15876 and /16678 
were measured. At the beginning of the outburst they devel- 
oped a P-Cygni type of profiles. From 2002 August absorp- 
tions affected the whole profile, mainly that of the /15876 line 
(Fig. 2). A narrow emission core without any wings and with 
absorption features at the top of the profile were observed. In 
addition, ratios of their fluxes differ significantly from those 
arising from net recombination. For example, we measured 
the ratio A = F(5876)/F(6678) ~ 1.2 H- 1.5, while flieo- 
retical inten sities from low-density plasma predict A ~ 3.6 
(e.g. Smits1 fl99^ . This difference has been ascribed to col- 
lisions and optical depth due to the metastability of the 2^S 
level, which affec t the He: level populations in a high den- 
sity plasma (e .g. lnegglll987HSchmidlll989at IProgaetal.l 
ll994clKin^don & Ferland 1119951) . Therefore we compared our 
flux ratios, F(6678) /F(5876)an d F(7065)/F(5876) with those 
modeled by Proga et al. I ill 994 for the case of high density 
nebulae in symbiotic stars (Appendix A, Fig. A.l). This sug- 
gests that the He -i- region is very compact with an average den- 
sity A^e Si 10'^ cm"-', radiating at ~ 20 000 K. Emissivity of 
such the plasma in the He i /15876 line corresponds to an effec- 
tive radius, 7?eif (He +) a; 5 Rq, which represents the radius of a 
sphere with constant density (see Appendix A in more detail). 

(v) Profile of the nebular [O iii] /15007 line showed asym- 
metrically extended core/wings with terminal velocities of 
about -300 and +250 km s"' . A sequence of profiles from 2000 
December 10, 1 1 and 12th, showed a creation of an extra emis- 
sion component on the violet side of the main core. In addition, 
weak emission components at about +380 km s"' accompanied 
the main core during the maximum. Their reality is supported 
by the S/N ratio of about 10 at these wavelengths and a rel- 
atively low level of the continuum from the giant (cf. Figs. 2 
and 7). If one associate these components to p olar blobs, as 
was interpreted for the classical Nova V 1 974 by Chochol et al. I 
{l997), the space velocit y of the eiecta would be larger than 
~1 600kms-i for / > 76° (HkQcalEoQ^- 

(vi) Evolution in the Raman scattered Ovi line at /16825 
provides information on the ionization structure around the ac- 
tive object. This emission was absent on our spectra from the 
maximum and after it. It appeared first on the 10/02/02 spec- 
trum as a very faint emission. During 2003 it became strong 
and broad (FW ZI 25+3 A) with 1+3 knots at the top of its 
profile (Fig. 4'). ISchmid1 il98 9b) suggested that this emission 
feature is due to Raman scattering of the O vi/il032 resonance 
line by atoms of neutral hydrogen. In Sect. 3.4 we investigate 
this process in detail. 

(vii) Transition to quiescence was characterized mainly by 
a significant increase in the hot object temperature (Sect. 3.3.3). 
As a result fluxes of highly ionized elements increased and new 
appeared (e.g. Ne vi lines. Fig. 5). At the end of 2003 a constant 
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Fig. 6. Dereddened FUSE spectra of ZAnd from the maximum of the optical brightness. It is composed from SiC2A 
(917-H1005 A), LiFl A (988^-1082 A) and LiF2A (1087H-1 181 A) channel spectra. Spectra from SiCl A (1080^-1090 A) and LiFlB 
(1080-^1188 A) have been used to complete the wavelength coverage. The solid sick line represents the Rayleigh attenuated 
(A^H - 7 10^' cm"^) blackbody radiation from the hot stellar source at 28 000 K (S ect. 3.3.2), while the dotted line is not attenu- 
ated. Most of the absorption features are due to H2 in the interstellar medium (cf. ICrowlevetar]l2004 . Note the P-Cygni type 
of profiles for C m, P v lines and the broad emission in the He 11 /11084 profile. 



brightness near to the values from quiescence and a very stable 
Ho- profiles (13/1 1/2003 and 18/12/2003) suggested the ZAnd 
active phase to be over. However, a strong wind from the hot 
star was still indicated by the P-Cygni profiles of the O vi line 
(Fig. 3), and in 2004 August we detected a new eruption in the 
LC (Fig. 1). 



drogen plasma. Then Eq. (1) takes the form dSkonallbOol 
Eq. 13) 

F, = F7"* (reff) + elnBMe-"''"'^ + ^ei(H, T,), (2) 



3.3. Spectral energy distribution 

Here we aim to reconstruct the SED prior to the maximum, 
at the maximum and during the post-outburst stage to under- 
stand better the observed photometric and spectroscopic evolu- 
tion throughout the outburst. 

3.3.1 . Model of the composite spectrum 

According to a simpl e photoionization model for symbiotic 
stars proposed by Seaauist et al. ] (Il984 (hereafter STB), the 
observed flux in the continuum, F{A), can be expressed as 

F(A)^F,(A) + Fi,(A) + FM (1) 

where Fn(/1) is the nebular component of radiation and Ff,(A) 
and FhiA) are stellar components from the hot star and 
the cool giant, respectively. We approximated radiation from 
the giant by a sy nthetic spectrum, F^'"'* (reff) (models from 
iHauschildt et al. i ri999). the hot stellar continuum by a black 
body at the temperature including the effect of Rayleigh 
scattering and the nebular radiation in the continuum by pro- 
cesses of recombination and thermal bremsstrahlung in the hy- 



where flj, - Id is the angular radius of the hot stellar source 
given by its effective radius R^^ and the distance d. Parameters 
determining the Rayleigh attenuation are the column density 
of H atoms, [cm^^], and the Rayleigh scattering cross- 
section for atomic hydrogen, cr^ [cm^]. Nebular radiation is 
in our approximation defined by the emission measure EM 
and the volume emission coefficient for hydrogen ei(H, Te) 
[ergcm-* s"' A"'], which depends on the electron temperature 
Te. The method of disentangling the composite spectra of sym- 
biotic stars was recently described bv lskopall 

To reconstruct the continuum of the hot object we used only 
the flux-points of the UBVR photometry and the far-UV contin- 
uum from ^1 188 to ^1 100 A. With the aid of the FUSE spectra 
we determined the column density of the atomic hydrogen by 
fitting the Rayleigh attenuated hot stellar radiation to the ob- 
served continuum blueward Ly/3 and Lya lines (Fig. 6). This 
reduced free parameters in the model to 9^, T^, and EM. As 
there was no ultraviolet spectrum available from the 2000-03 
outburst, determination of Tp is more uncertain. To satisfy the 
condition for an appropriate fit ( Skopal 2005, Sect. 3.3) we es- 
timated the uncertainty in Th within the range of 500-H 1 000 K, 
while Te can run between about 30000 and 40000 K. 
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3.3.2. SED at the plateau stage and the maximum 

To reconstruct the SED during the time of the last plateau in 
the LC we used the optical spectrum from 06/11/00 and the 
FUSE observation from 16/1 1/00. For the maximum, the cor- 
responding data are from 11/12/00 and 15/12/00, respectively. 
A comparison of our solutions with observations is plotted in 
Fig. 7 and the corresponding parameters are in Table 2. 

The change in the SED from the plateau at C/ ~ 9.0 to the 
maximum {U ~ 8.3) was due to an increase in both the hot 
stellar and the nebular emission. These components of radia- 
tion contributed more to the U band than to the BVR region. 
Therefore we observed a larger increase of the star's brightness 
in U (AU ~ 0.7, AB ~ 0.4, AV ~ 0.3, AR ~ 0.3 mag), i.e. the 
U-B index became more blue. The temperature of the hot stel- 
lar source, Th, increased from about 26 000 K to about 28 000 K 
at a constant effective radius (Table 2). Accordingly, its lu- 
minosity increased from 3 800 to 5 300 {d/1.5 kpc)^ Lq during 
this transition. However, the emission measure also increased; 
significantly by a factor of 2.2. 

The observed hot stellar radiation is not capable of pro- 
ducing the nebular emission (too low T^). This implies that 
the hot ionizing source cannot be seen directly by the outer 
observer, but can be 'seen' by the surrounding nebula, which 
converts its radiation into the nebular emission. This suggests 
a disk-like structure of the optically thick material, whose 
outer rim occults the central hot star in th e dkection of the 
observer (also Sect. 5.3.5 of ISkonallEoOsI) . For comparison 
we show the SED from the previous 1986 optical maximum, 
which is of the same type. In such the case the temperature 
of the ionizing source, T'^^\ can be determined only indirectly 
from the nebular radiation , which m i rrors i ts radiative prop- 
erties. For example, .Sokoloski et al. ] ll2005 l) determined T/^ 
= 95 000 K from fluxes of Hj8 and Hen 4686 emission lines 
observed at the maximum. However, scaling the hot radiation 
from the unseen ionizing source to the observed, significantly 
cooler, pseudophotosphere yields an erroneous (too high) lu- 
minosity of the white dwarf. According to our SED model 
the total luminosity produced at/around the white dwarf dur- 
ing the maximum is Lt S ih + = 6 900 (dl 1 .5 kpc)^ Lct) 
(Table 2 here and Sect. 5.3.6 of Skopal '200?, in more detail), 
but Sokol oski et al, ( 2005) derived Lh = 9 800 (t// 1 .2 kpc)^ Lq 
for £b-v = 0.27 (= 21 500 (c//1.5 kpc)^ L© for £b-v = 0.30) 
by scaling the 95 000 K radiation of the ionizing source to the 
observed FUV fluxes from the warm pseudophotosphere. Also 
the best model SED of the FUV and optical continuum with the 
fixed stellar component radiating at T^^^- - 95 000 K is not sat- 
isfactory (x^56.2). For comparison, our models in mid panels 
of Fig. 7 have and;if=13.6, respectively. 

3.3.3. SED at the post-outburst stage 

Here we selected the last observations made by the FUSE 
on 04/08/03 and corresponding UBVR photometric measure- 
ments. The slope of the B, V and R fluxes has an opposite trend 
to that from the maximum - the fluxes follow rather the SED 
of the giant, which dominates the spectrum from the B band 
to longer wavelengths (Fig. 7). The U flux-point does not obey 
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Fig. 7. The UV/optical/IR SED during active phases of Z And. 
Top: At the 1985-maximum. Middle two panels: Prior to and 
at the maximum of the 2000 outburst (positions in the LC are 
shown in Fig. 1). Bottom: At the post-outburst stage in the sum- 
mer of 2003. The resulting spectrum (solid thick lines) is given 
by superposition of nebular component(s) (dashed lines) and 
two stellar components, from the hot stellar source (solid thin 
Unes) and the giant (dotted lines). IR fluxes (O) are from quies- 
cence and far-UV fluxes (o, +, Sect. 2.2) are near-simultaneous 
to the optical observations. 



this trend. It is higher, the ratio Fu/F^ ~ 1.78 (i.e. U - B - 
-1.16). This reflects a dominant contribution from the nebula 
in the near-UV/U. As a result the radiation from the hot star in 
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Table 2. Fitting and derived parameters from the SED during the maximum of the 1985 and 2000 outburst. 



Date 






pert 




log(nH) 




EM 


Ln 


F(/5£- Optical 


[K] 






[Lq] 


[cm-^] 


[K] 


[cm-3] 


[Lq\ 


24/12/85 


19000 


3.4 10-" 


2.3 


590 


23.34 


18 000/40000 


2.1 10*" 


425 


16/11/00-06/11/00 


26000 


4.6 10-" 


3.1 


3 800 


21.85 


37 000 


3.9 10''° 


740 


27/11/00-11/12/00 


28 000 


4.7 10-" 


3.1 


5 300 


21.85 


38 000 


8.8 10*" 


1600 


04/08/03-31/07/03 


120000 


1.510-'2 


0.10 


1800 




25 000 


1.3 10«" 


280 




160000 


1.3 10-'^ 


0.086 


4400 











the optical had to be small relatively to other sources here. This 
suggests a high temperature, Th ^ 100 000 K. We adopted Tj, = 
120 000 K to satisfy basic p/iofo-ionization model of symbiotic 
stars (STB). At this temperature the hot star luminosity is just 
capable of producing the observed EM. From this point of view 
Tb > 120 000 K. Finally, we complemented contributions from 
the giant and the hot stellar source by a nebular contribution 
with Te = 25 000 +5 000/- 7 000 K and EM = 1.3 lO'-^cm-^ 
Large uncertainties in are given by the absence of ultravio- 
let data. Model parameters are in Table 2. 

[Sokoloski et al. (2005) determined = (160 000 + 
35 000) K from emission line fluxes and Rf = (0.06 + 0.01) x 
1.47 = (0.088 + 0.015) /?0, Lh = (2000 + 500) x 2.20 = 
(4400 ± 1 100) Lq from the FUV continuum fluxes. The factors 
1.47 and 2.20 convert their original quantities calculated for 
Eb-y - 0.27 and d - 1.2kpc to those corresponding to /iB-v = 
0.30 and d-\.5 kpc we adopted in this paper. In this case phys- 
ical parameters derived from the SED (Table 2) are the same as 
those given by optical line fluxes in combination with the FUV 
continuum fluxes. This reflects an optically thin regime in 2003 
- in contrast to the situation from the optical maximum. 

3.3.4. SED as an indicator of tine mass outflow 

The total emission measure of the nebular continuum at the 
maximum, EM ~ 6 10''" cm-^ (Table 2), was a factor of 
5-^10 larger th an in quiescence (iNussbaumer & Vogel I !l989: 
ISkopallf2005' ) and a factor of 3 larger than at the m aximum 
of the 1985 eruption (Fernand ez-Castro et al. lll995l) . During 
quiescent phases the nebula originates mostly from ioniza- 
tion of the giant's wind. It is characterized by a large param- 
eter X > 10 in the STB model (e.g. Fernandez-Castro et al. 
ll988HSkonal1l200ll) . which corresponds to a very open neb- 
ula (i.e. in major part a particle bounded). Therefore the sig- 
nificant increase of the nebular emission during the maxima 
of the outbursts can result from a supplement of new emit- 
ters into the particle bounded nebula, accompanied probably 
by an increase in the flux of ionizing photons from the hot 
source. The former can be caused by an increase in the mass- 
loss rate from the hot object. The broad Ha wings (Fig. 4, 
Sect. 3.5) suggest t he ma ss-outflow to be in the form of a stel- 
lar wind. ISkot)all 1I2OO6 ') estimated the corresponding mass- 
loss rate at the maximum as ~ 1 -H 2 10'^ Mq yr^\ which is 

a factor of ~ 3 larger than that from the giant (see Table 3 in 

I ■ ° ° 

ISkop al 2005). The presence of a rather high-temperature neb- 
ula (T^e = 30000 -H 40 000 K, Table 2) requires a large contri- 
bution from collisions (.Gurzadvan ■■1997.) . which is consistent 



with the strong mass-outflow. Also lSokoloski et al. I ( l2005'l) as- 
cribed the increased X-ray emission at the optical maximum to 
the shock-heated plasma as a consequence of the mass ejection 
from the white dwarf into the dense symbiotic nebula. Also 
the transient increase of the nebular emission during the low- 
amplitude 2002/03 rebrightening (bottom panel of Fig. 7) was 
also connected with a release of new emitters into the particle 
bounded nebula (see Sect. 3.6). Contrarily, during the plateau 
stages the SED did not change considerably as suggested by 
the stability of colour indices (Sect. 3.1). This constrains the 
hot object to be stable in the luminosity, temperature and also 
the mass outflow, because more emitters would give rise a sur- 
plus of the emission measure and thus the star's brightness. 



3.4. Direct and Raman scattered O \/\A1032 line 



The Raman scattering process a nd its possible appl ications in 
astrophysics was introduced bv INussbaumer et al. I ([l989). In 
this process a photon excites an atom from its ground state to 
an intermediate state being immediately stabilized by a transi- 
tion to a true bound state by emitting a photon of a different fre- 
quency. "Schiiiid" (798%) identified the .16830 and /17088 emis- 
sion bands observed in many symbiotic stars as to be due to 
Raman scattering of the O vi resonance doublet /11032, /11038 
by neutral hydrogen. Later on this process has been investi- 
gated by many authors (e.g..Schmid 1996: Harries & HowarthI 
ll 994ISchmid et al. 111999^ iBirriel et al. Il200nl) . The efficiencv 
of this scattering process and profiles of both the direct and 
the Raman-scattered O vi lines together with polarimetric mea- 
surements provide i nformation on the scattering geometry in 
symbiotic systems (Schmid 1998). Previous studies referred 
exclusively to quiescent phases of symbiotic stars. For Z And, 
iBirriel etal.i ri998) revealed firstly a strong evidence for the 
Raman effect and determined its efficie ncy to ~ 7 % for /1 1032 
on the basis of the HUT observations. Schmid et al. ( 1999) de- 
rived the efficiency of this conversion to 4^6% and found that 
the Raman /16825 profile in Z And is by 30% broader than the 
original Ovi/11032 line on the violet side of the profile. They 
used observations from the ORE E US mission. 



Here we analyze our series of Raman /16830 lines ob- 
tained along the major ZAnd outburst (Fig. 4) together with 
the archival FUSE spectra covering the same period. 
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Table 3. Dereddened fluxes for the Ovi/}1032 line and 
its Raman scattered counterpart at /16825 in units of 
10"'^ erg cm"- s"'. The efficiency, rj, was derived according to 



flie relation (3). 


Date 


Julian date 
24... 


Phase 


^Ovi 




'7 

[%] 


16/11/00 


51865.0 


0.16 


50.6 






27/11/00 


51876.0 


0.18 


44.6 






10/12/00 


518894 


0.19 


45.4 i 


<2t 


< 29 


11/12/00 


51890.4 


0.19 


45.5 i 


<2t 


< 29 


12/12/00 


51891.3 


0.20 


45.6i 


<2t 


< 29 


15/12/00 


51894.0 


0.20 


45.8 






05/01/01 


51914.9 


0.23 


74.5 i 


<2t 


< 18 


27/05/01 


52057.3 


0.42 


270. i 


<2t 


< 5 


22/07/01 


52113.3 


0.49 


347. 






30/09/01 


52183.3 


0.58 


308. 






10/02/02 


52315.9 


0.76 


972. i 


4.31 


2.9 


05/07/02 


52461.5 


0.95 


1700. 


8.93 i 


3.5 


07/08/02 


52494.3 


0.99 


1520. i 


9.97 


4.3 


22/10/02 


52570.5 


0.09 


1103. 


16.2 i 


9.7 


02/02/03 


52673.5 


0.23 


1303. i 


24.7 


12.5 


31/07/03 


52852.5 


0.47 


1650. i 


13.7 


5.5 


04/08/03 


52856.8 


0.47 


1658. 


13.8i 


5.5 


02/09/03 


52885.5 


0.51 


1714. e 


14.4 


5.6 



i/e ~ interpolated/extrapolated value, 
t - no detection, upper limits only 



3.4.1. Raman scattering efficiency 

The scattering efficiency for the Raman process is defined as 
the photon ratio, 77(0 vi) - A^Ram/A^ovi, between the Raman 
scattered and the initial O vi line component, i.e. the fraction 
of emitted O vi photons converted to Raman photons. As we 
measured only the Raman /16825 line we investigate here the 
conversion /1 1032 — > /16825. In the simplest case the efliciency 
of this transition can be written as 

'?(Ovi) = 6.614^, (3) 

■^Ovi 

where the factor 6.614 = il031.92//16825.44. In the real case 
77 depends on the geometry of H i atoms within the binary and 
the binary position with respect to the observer ( Schmid.,1996,) . 
Thus the relation (3) represents just an empirical measure of 
the line photon ratio we call efliciency. Table 3 and Fig. 8 sum- 
marize the results. As the FUSE and our optical observations 
were not simultaneous, we interpolated FUSE fluxes to dates 
of optical observations and vice versa. Below we comment on 
the variation in 77 along the outburst as follows: 

(i) At the initial stages (2000-01) the not-measurable 
Raman's band could be a consequence of a very low O vi /il032 
flux. Assuming the conversion efficiency of a few percent, the 
incident flux of 45x 10"'^ erg cm"^ s"' (Table 3) produces a few 
times lO"'"* erg cm"^ s"' at /16825, which implies the height of 
the Raman feature above the local continuum only of a few 
times 10 '"^ergcm^^s which is under the detection. On the 
other hand, the signal to noise ratio S/N~20 around the Raman 
line suggests a limiting flux of ~ 1 x lO"'-' ergcm"- s"' A"' 
to recognize a real detail in the spectrum. Then, assuming the 



Raman line to be flat and broad as ~20 A, its flux we cannot 
detect is $ 2x10 '^ erg cm"^ s ' . Corresponding upper limits 
for efficiencies 77 are relatively large, because of low initial O vi 
fluxes at the optical maximum (Table 3, Eq. (3), Appendix B). 

(ii) In the period of 2002, prior to the spectroscopic con- 
junction, the emergence of the Raman features and the en- 
hancement of the O VI fluxes (Table 3, Fig. 8) require a relevant 
extension of the scattering region to get the measured efficien- 
cies of 3-h4%. Note that the flux increase of the /11032 line 
itself would result in a decrease of the conversion efficiency 
(Eq. (3)). A comparison of the d erived efficienci es to a grid of 
model calculations published bv'SchmidI (Il996l) suggests that 
a kind of his XB3 model could be relevant, i.e. nearly one third 
of the sky, as seen from the O vi ionized zone, is covered by the 
neutral scattering region. 

(iii) A significant increase in the Raman conversion effi- 
ciency during the first quadrature after the spectroscopic con- 
junction (2002 autumn - 2003 spring, orbital phases 0-^0.25, 
Fig. 8) reflects an enlargement of the neutral region seen from 
the Ovi zone. The observed Ovi/11032 flux decreased, while 
its Raman scattered counterpart at /16825 increased (Table 3). 
Here the efficiency of 10-h12% are similar to that of the XCi 
model of Schmid's (1996) simulations, in which at least one 
half of the 'O vi sky' is occupied by atoms of neutral hydrogen. 
This change in distribution of the neutral material in the binary 
could be connected with an enhancement of the hot star wind 
near the equatorial plane as a result of dilution of the optically 
thick disk-like shell during this period (Sect. 4). According to 
hydrodynamical calculations that include effects of the orbital 
motion, a dense S-shaped wind-wind collision region can be 
formed between the stars where the fast wind pushes the giant 
wind (e.g. Fig. 3 of lWalderlll995l) . This interaction zone can 
divide the binary space into two approximately equal parts of 
the neutral and ionized wind. Such a distribution of neutral gi- 
ant's wind in the active binary satisfies conditions required by 
a larger efficiency of the Raman scattering process. 

(vi) During the post-outburst period (2003, orbital phase 
~0.5), observed Ovi fluxes increased, while the /11032 — » 
/i6825 conversion efficiency decreased. This indicates a de- 
crease of the size of the scattering region as 'seen' by the O vi 
photons. De rived efficiencies of ~ 5 % are similar to those from 
quiescence llBirriel et al. ll998HSchmid et al. Il999l) . 

3.4.2. Comparison of line profiles 

Evolution in the Raman scattered Ovi /1 1032 line along the 
outburst is shown in Fig. 4. A comparison of both direct and 
Raman scattered line is plotted in Fig. 9 for near-simultaneous 
optical/far-UV observations. Table 4 summarizes their basic 
parameters {RVc, FWHM). Well-measurable Raman /16825 
line developed from QllQ^IQl on our spectra, i.e. nearly two 
years after the beginning of the outburst. Profiles consist of 
a central peak at -6832 A and a violet and red shoulder at 
~ 6823 A and ~ 6837 A, respectively. On 02/02/03 the red adja- 
cent component developed in a separate peak of a comparable 
intensity to the main peak. The directly measured Ovi /1 1032 
line was characterized by a strong central emission shifted 
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Orbital phase 
0.25 0.5 0.75 1 1.25 1.5 



Table 4. Line center (RVc) and line width (FWHM) for direct 
and Raman scattered O vi /11032 line in km s ' . 
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Fig. 8. Evolution of the Raman conversion efficiency /11032^ 
/16825 along the outburst of Z And (full diamonds). Open trian- 
gles connected with dotted line denote upper limits for too low 
Raman fluxes that were beyond our detection (see the text). 
Data are from Table 3. 



by about -1-58 km s"' with respect to the systemic velocity. 
Such a l arge redshift was mea sured also during a quiescent 
phase by ISchmid et al. I (Il999l) who explained this effect by 
radiative transfer in strong resonance lines throughout an ex- 
panding medium. During the active phase, this line absorp- 
tion/scattering effect attenuated the violet wing from about -1-50 
to -350 km s ' and produced a P-Cygni type of the profile 
(Fig. 3). This suggests its origin in the densest part of the hot 
star wind. Other variations along the active phase (Fig. 9) are 
discussed below. 

We compare the line profiles in the radial velocity (RV) 
space considering the principle of energy conservation of the 
Raman scattering process. According to iNussbaumer et al.1 
( Il989l) a broadening AAi of the incident line Ai is transferred 
to a broadening AAf of the scattered fine Af as 



AAf 
IT 



~A~J' 



(4) 



In the RV-scale and for the Ovi(/I1032 /16825) scattering 
process, Eq. (4) can be written as ARVoyi = A/?yRain/6.614. 
Further we assume: (a) The O vi emission region is connected 
with the wind from the hot star, (b) The scattering region is rep- 
resented by neutral atoms of hydrogen of the wind from the gi- 
ant, (c) Between the binary components there is an interaction 
zone of the two winds. It separates a low density (filled with 
material from the hot star) and a high density region (filled with 
material from the gi ant). The zone c an be spiral-shaped due the 
orbital motion (e.g. IWalderlfl995l) . With the aid of a general 
discussion on the O vi line profiles of 6 objects as given by 
ISchmidet al"l(ll99 9V we propose the following explanation of 
our observations (Fig. 9): 

(i) The first measured Raman /16825 line (10/02/02, (p = 
0.76) was compared with the direct O vi line from 30/09/01 ((p 
= 0.58). Unusually, with respect to later observations, both pro- 
files nearly overlapped each other. So there were no significant 
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motions of the scattering region relative to the direction of the 
incoming O vi photons. This implies that the O vi photons could 
not impact the neutral wind at the orbital plane moving to- 
ward them. Otherwise the ~ -40 km s ' shift of t he Raman line 
relative to the original line had to be observed ( Schmid et al. I 
■1999.) . According to our interpretation of the hot object struc- 
ture (Sect. 4), a disk-like material blocked the O vi photons at 
the orbital plane, which thus could 'see' only outer parts of the 
giant's wind giving rise to only a small redshift in the /16825 
profile. Note that the major part of the O vi ffux is created in 
the densest regions of the wind encompassing the hot star (cf. 
Appendix B). Another peculiarity is that the violet wing of the 
/1 1032 line was not so strongly absorbed in the wind medium 
as we observed later during the outburst. As a result the total 
profile was shifted only to — 1-35 km s"' (Table 4). A shallow 
minimum in the /1 1032 line at ~ -)-18kms"' could be caused 
by a self-absorption near to the hot object. Note that the ra- 
dial velocity of the hot component orbital motion at the time of 
the observ ation ((p = 0.76) is comparable t o this shift (param- 
eters from Mik olaiewska & Kenvo n 1996), which thus places 
its origin at/around the hot star, in agreement with our assump- 
tion (a) above. 

(ii) Next available pair of Ovi/11032 and /16825 lines was 
taken just prior to the inferior conjunction of the red giant {tp 
= 0.95 and 0.99). The fine ffux at /11032 was the largest from 
our sample and its profile was characterized by a narrow red- 
shifted core (RVc ~ -1-67 km s"'), while its scattered counter- 
part at /16825 was rather faint. The profile of the initial Ovi 
line was probably modified by an occultation of the O vi neb- 
ula, which thus reduced its violet wing for the outer observer. 
Note that from p ~ 0.94 a decrease of th e optical light du e to 
the eclipse effect was observed (Fig. 2 of Skopal 2003a). On 
the other hand, the O vi zone was not subject to occultation for 
the scattering region spread within the H i zone in the binary. 
So this could 'see' the un-afffected Ovi/11032 line profile and 
thus convert it to a standard type of the /16825 line profile as 
measured at other orbital phases later on. The significant in- 
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Fig. 9. Comparison of the Ovi/11032 line profiles with their 
Raman scattered counterparts along the outburst. Compared is 
also the [Fevn]/i6087 line profile shifted by -1-45 km s"' (see 
Sect. 3.4.3). The O vi line flux was scaled to that of the Raman 
line. 



from the giant moving against it. This lead to the appearance 
of the violet adjacent component in the /16830 profile and an 
enhancement of the 77(0 vi) efficiency. 

(iii) After the spectroscopic conjunction (the far-UV and 
optical observation made at 1^ = 0.09 and 0.23, respectively; 
Table 3), during the 2002/03 transient increase in the star's 
brightness, we observed the largest efficiency of the /1 1032 — » 
/16825 transition, well developed three-peaked structure of the 
Raman line and a relatively broad line core of the direct O vi 
line with superposition of a wind component at ~ -1- 180 km s ' 
(Fig. 9). The Raman line was broader by the blue adjacent com- 
ponent. Its wing was shifted by about ^Okms"' relatively to 
that of the /11032 line. This emission is produced by the Hi 
atoms of the giant's wind, which are located along the binary 
axis and move toward the O vi region. The absorption compo- 
nent in the emission core of the Ha profile at the position of 
X -45kms"' (Sect. 3.2 (ii)) supports this interpretation. In ad- 
dition, a dense neutral region could be created as an enhanced 
hot star wind pushes the neutral giant's wind between the bi- 
nary components. This can increase the efi'ective number of the 
hydrogen atoms and thus also the efficiency 77(0 vi). The outer 
part of the neutral wind and the Hi atoms in the collisional 
zone moving from the incoming O vi photons then extend the 
red side of the Raman line profile. The red wind component in 
the /1 1032 profile has no Raman counterpart. This portion of 
the O VI emission has an origin at outer part of the hot wind, 
i.e. further away from the scattering region, which reduces the 
scattering efficiency. 

(iv) Similar profiles were observed also later at the orbital 
phase if - 0.5 (nearly simultaneous observations, bottom panel 
of Fig. 9). This suggests the difference between the direct and 
the Raman line profile to be of the same nature as we discussed 
above in the point 3. 

The shift between the red wings of the Raman-scattered 
and the direct Ovi line during a quiescent phase was larger 
than we observed during this activity. FWHM of the Raman 
line from the activity exceeded this quantity from quiescence 
by about 10H-30kms whil e the width of the incident Ovi 
line d id not change so much ( Schmid et al. ll999tlBirriel et al. I 
l2000t Table 4 here). This reflects a different kinematics of the 
scattering region during active and quiescent phases. Because 
of a strong wind from the hot star developed during the ac- 
tive phase, this could be due to the wind-wind interaction zone 
in the binary. Therefore, in the following section, we discuss 
possible connections between profiles of the [Fe vn] /16087 and 
Raman /16825 line. 

3.4.3. Raman ^6825 and [Fe vii] /i6087 profiles 



crease of the direct O vi line flux at approximately unchanged 
or slightly increased conversion efficiency constrains also an 
enhancement of the H i atoms taking part in the scattering pro- 
cess (Eq. 3). The former resulted from the temperature increase 
of the ionizing source (Sects. 3.2 (vii) and 3.3), while the latter 
could result from a dilution of the cooler disk-like pseudopho- 
tosphere. As a result the wind from the hot star could appear 
close to the orbital plane and thus met a massive neutral wind 



A comparison of [Fe vii] /16087 and Raman /16825 line profiles 
is shown in Fig. 9. The former was shifted by +45kms"' to 
eliminate the redshift in the Raman line (Sect. 3.4.2). Similarity 
of both profiles is obvious, mainly in 2002 and 2003 when 
also their line fluxes were comparable. Their profiles - a large 
FWHM of x!l20-140kms"', steep wings with a multi-peak 
structure at the top - suggest them to be due to kinematics of 
emitting regions. Wings of both the lines practically overlap 
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each other at least up to the half of the maximum. Striking is 
the nearly equality of the extended red wings with a pattern of 
the red adjacent component (bottom panels of Fig. 9). 

These common signatures of both line profiles suggest a 
similar kinematics for the responsible emitting regions, but dif- 
fering drastically in conditions of ionization - one region has 
to contain Hi atoms at high density, while the other one the 
Fe^'' ions at a low density. This can be satisfied if these lines 
are formed within the wind-wind interaction zone. The Raman 
line can be created at/around the interaction surface from the 
giant's side, where its neutral wind is pushed by the hot star 
wind, while the [Fevn] emission can arise in a collisionally 
ionized boundary layer between the two winds. The observed 
change in the line width at ^ ~ 0.25 (02/02/03, a maximum 
width) and ~0.5 (31/07/03, a narrower width) is consistent 
with the main wind flow velocity aligned approximately per- 
pendicular to the line connecting the stars. Previous works also 
support such origin of [Fe vii] lines. Wallerstein et al. ' ('1984') 
showed that line profiles of forbidden lines from highly ionized 
atoms could result from interaction nebula of wind-wind colli- 
sion in a binary system. The observed profile of the [Fe vii] line 
in V1016Cyg was very similar to our observations from 2003 
(see their Fig. 4). Also hydrody namical calculations of the col- 
liding zone of the two winds by Nussbau mer & Walder 1 lll993l) 
showed a wide variety of resulting line profiles, some of which 
are similar to our observations (e.g. the overall broadening with 
steep sides of profiles). 

Accordingly the observed broadening of the Raman line 
profile together with profiles of highly ionized forbidden lines 
can provide a diagnostic tool in studying motions in the wind- 
wind interaction zone during active phases. Note that a Raman 
photon may escape easily from the neutral region, because this 
elastic process cannot be repeated (in contrast to Rayleigh scat- 
tering) and thus reflects motions in the densest (a small cross- 
section) outer layers illuminated by O vi photons. 

3.5. On the origin of tine broad Ha wings 

iNussbaumer et al.1 (Il989l) pointed out that the broad Ha 
wings in symbiotic stars could be formed through Raman 
scat tering of Lyyg photon s by atomic hydrogen. Lee (2000) 
and iLee & Hvung elaborated a relevant model. They 

showed that the wing profile formed through the Raman 
scattering can be approximated by the curve proportional to 
f (Av) PC Av~^ in the velocity co-ordinates. On the other hand, 
LSkopal et al. | (l2002all derived an expression of the same type 
(see their Eq. 1 1) to fit the broad components of the hydrogen 
line profiles that form in an optically thin and spherically sym- 
metric stellar wind. Therefore it is not possible to distinguish 
the wind emission from the scattered features in the Ha wing 
directly. In addition, the observed LyyS profile can be influenced 
by other processes in the circumstellar medium. 

In our case the observed Ly/3 profile is not able to produce 
the Ha wings through the Raman scattering process. The width 
of the scattered Ha wings is proportional to the width of the 
LyyS line as AvLy/j = AvHff/6.4. Then the full extension of Ha 
wings ^ 4 000km s"' (Fig. 4) requires AvLy/j > 625 km s"' 



or A/lLy/s > 2.1 A, which is a factor of ~2 larger than the ob- 
served FWZl(Lyy6) < 1 .0 A. In the case of an optically thick 
Rayleigh scattering region the Ly/3 photons may undergo many 
Rayleigh scatterings until they are converted into Raman pho- 
tons and escape. Thus in such the special case only Raman fea- 
tures can be detected, but no original Ly/? emission. The real 
case is complicated with the presence of a strong and variable 
geocoronal component of LyjS, which is difficult to isolate from 
the circumstellar one. As a result we have to investigate other 
possibilities to identify the main source of radiation co n tribut- 
ing to the broad Ha wings. For example, llkeda et al 1 ll2004 
found that polarization profile of Ha does not agree with that 
of the Raman /16830 line for Z And on 25/10/02 (i.e. during the 
2002/03 brightening). This finding suggests that the Ha line 
does not include the Raman-scattered component that origi- 
nates from Ly/3. Also a relationship between the emission from 
Ha wings and the emission measure of the symbiotic nebula 
supports the origin of the broad Ha wings in the ionized stellar 
wind (Skopal 2006). 

Therefore we suggest that the extended Ha wings, which 
developed during the investigated active phase of ZAnd, are 
due to kinematics of the ionized hydrogen rather than the 
Raman scattered Ly/3 photons. However, this problem deserves 
further investigation. 

3.6. Optical rebrightening in 2002/03 

After 2002 September/October, when the hot object was aris- 
ing up from the eclipse, a rebrightening in the LC (AC/ ~ 0.6, 
AB ~ 0.4, Ay ~ 0.3) was observed for about 5 months. Our es- 
timate of its magnitude increase takes into account the depth of 
the eclipse and the brightness level prior to it ( Skopal 20 03al) . 
During this event we measured an additional ex tension of the 
Ha wings on our spectrum from 02/02/03 and iTomov et al. I 
(2005) observed such broadening also in the emission wings of 
the He ii /14686 line profile with FWZI~2400km s ' . 

The observed indices from the maximum (around 
15/12/2002) were (U- B)™" = - 0.70 and (B - V)™" = 0.85. If 
we deredden these observations with E(B - V) = 0.30, correct 
for emission lines (we adopted At// = 0.25, AB/ - 0.45 and 
AVi = 0.18 coii'espondi ng to a st andard emission line spectrum 
for classical symbiotics dSkopal 1120 03b)') and subtract contribu- 
tions from the giant, Fg(B) = 1.3 x lO"'-' erg cm"^ s"' A"' and 
Fg(y) = 2.9 X lO^'^ergcm^^s"' A"' (see Fig. 7), we get in- 
dices of the extra light that causes the optical rebrightening as 
(U - B)™^,= - 1.27 and (B - V0™ra=0-59. These quantities re- 
flect a large contributi on from the nebula rad iating at high elec- 
tron temperature (e.g. iMenzies et aLlll982[) . We demonstrate 
this case in Fig. 7. The UBVR flux-points from the maximum 
of the 2002/03 rebrightening suggest the same type of the SED 
as in the mid-2003, but with a factor of ~2 stronger contribu- 
tion from the nebula (^^^(rebrightening) a; 2.6 x 10^" cm"^). 

Assuming a particle-bounded nebula (i.e. geometrically 
open in the sense of the STB model), a maximum increase 
of its radiation at a constant flux of ionizing photons can be 
produced by an injection of new particles into it, that will con- 
vert a relevant amount of the ionizing photons into the nebu- 
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lar radiation to yield a radiation-bounded (closed) nebula. The 
particle-bounded nebula in Z And during 2002-03 is supported 
by a high quantity of the hydrogen ionizing photons of ^ 1 .2 x 
10^^ s"' given by the hot object luminosity (Lh > 1 800 L0) 
and its temperature (Th > 120 000 K, Table 2). A large value 
of the parameter X (x; 14) during quiescence corresponding to 
an open nebula was noted by Fernandez-Castro et al. ( 1988). 
Within this scenario we interpret the small optical rebrighten- 
ing as a consequence of a dilution of the optically thick shell 
around the hot star This event supplied additional emitters into 
the particle-bounded nebula, which can convert a relevant part 
of the far-UV radiation through f - b and f - f transitions to 
the optical. By other words, the material became more trans- 
parent and thus was subject to ionization by the hot central star 
even at the orbital plane, which led to a transient increase of 
the nebular emission. We add that similar brightening was also 
observed in the LC of AX Per at its tr ansition from the ac tive 
to quiescent phase (Figs. 3 and 10 of ISkopal et al. lEoOlh . A 
dilution of an optically thick shell in this case was directly con- 
firmed by the change of the naiTow eclipses into very broad 
minima occurring just after the event (cf. Fig. 3 there). 

ISok oloski et al. I ( Eo05 ) interpreted the 2002/03 rebright- 
ening as a result of a decrease in (from around 150 000 K 
to approximately 120 000 K) due to a slight expansion of 
the white-dwarf photosphere. They based this view on that 
both the FUSE and VLA fluxes were unusually low during 
the optical rebrightening and the U - B index varied from 
about -0.45 to -0.8 from the beginning to the end of the 
event. We put following arguments against this interpretation: 
(i) The decrease in the FUV fluxes and Th implies a de- 
crease in both the stellar and the nebular component of radi- 
ation from the hot object. Note that such the small temper- 
ature decrease at Th > 100 000 K induces a negligible (i.e. 
not m easurable) change in the SED (e.g. Munari & Buson 
Il992h . So, we should observe a decrease of the optical bright- 
ness, (ii) There is an additional attenuation of the UV con- 
tinuum that is well pronounced for eclipsing systems at po- 
sitions around the inferior conjuncti on of the giant. Thi s ex- 
tinction process was first noticed by Miirset et al. '^('1991') and 
[pumm et al. ( 1999) analyzed it for RWHya. Skopal (20051 
found this effect to be present in aU lUE spectra of ZAnd in 
the range of phases tp ^ + 0.15. Therefore the 'unusu- 
ally' low FUSE fluxes from 22/10/2002 ((p = 0.093) result 
from this eff'ect. (iii) The observed fading in the radio fluxes is 
consistent with that in EM: £'M(maximum)/£'M(2002/03) ~ 
^^5GHz(inaximum)/F5GHz(2002/03) within th e corresponding 
uncer tainties (Table 2 here and Table 5 of ISokoloski et al. I 
l2005h . (iv) A redder U - B index at the beginning of the re- 
brightening (-0.45) and a more blue values at/after the maxi- 
mum (-0.8) could qualitatively reflect variation in optical prop- 
erties of the ionized medium during/after the shell dilution. 



4. Structure of the hot object during the outburst 

4.1. The stage at/around the optical maximum 
This stage was characterized by: 
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(i) The two-temperature-type of the hot object spectrum. 
The cooler one was produced by a relatively warm stellar 
source radiating at 20 000^30 000 K (Fig. 7, Table 2) and the 
hotter one was represented by the highly ionized emission lines 
(e.g. Hen, Ovi, [Fevii]) and a strong nebular continuum ra- 
diation (EM ^ 10^° cm"^. Table 2). The former is not capa- 
ble of producing the observed amount of the nebular emission 
and thus the latter signals the presence of a hot ionizing source 
(k, 10^ K) in the system, which is not seen directly by the outer 
observer 

(ii) Signatures of a mass-outflow at moderate veloci- 
ties (~ 100-H200kms"'), as indicated by the P-Cygni profiles 
(Sect. 3.2 (i). Figs. 2 and 6; also Sokoloski et al. 2005), and at 
very high velocities (^ 1000-^2000 km s"') having a pattern in 
broad emission wings of Ha (Fig. 4). 

The two-temperature spectrum and the two-velocity type 
of the mass-outflow from the hot object suggest an optically 
thick disk-like structured material surrounding the central hot 
star and expanding at moderate velocities at the orbital plane, 
while at higher latitudes a fast optically thin wind escapes 
the star. Due to a high orbital inclination the outer observer 
can see just the optically thick matter of the disk-like shaped 
warm pseudophotosphere, whose outer flared rim occults the 
central ionizing source. Only the suiTounding medium from 
above/below the disk can 'see' directly the ionizing source and 
thus convert its radiation to the nebular emission. As a result, 
the temperature of the ionizing source, T ^'^-, ca n be obtained 
only through emission lines ( Sokolos ki et al. ]|2005. derived 
T^^- ~ 95 000 K), while the observed UV/optical continuum 
determines the temperature of the stellar source (i.e. the warm 
pseudophotosphere), ~ 27 OOOK (Sect. 3.3). Thus this stage 
was characterized by 

« n^ (5) 

Also the low fluxes of highly ionized elements, whose profiles 
were shifted blueward (Sect. 3.2 (iii). Fig. 2 for Hen, Fig. 4 
for [Fevii]), are consistent with the disk-like structure of the 
hot object. A small size of the corresponding ionized zone (see 
Appendix B) produced relatively weak emission and the op- 
tically thick disk-like extended source, inclined at 75 - 80° 
J ' 11 I 

(|Skopal 2003a), blocked a fraction of radiation from an out- 
flow at its back side in the direction of the observer. As a result 
the blue-ward shifted emission dominated the profile. This ba- 
sic structure of the active object persisted probably until about 
2002 August, when the eclipse effect in the LC was the last 
indication of the pre sence of the op tically thick disk-like shell 
around the hot star ( Skopal 2003a). 

Finally, we note that such the structure of the hot object 
is common for all active S -type symbiotics with a high orbital 
inclination dSkopal l2005l Fig. 27). 

4.2. Evolution after 2002 August 

During this period we indicated a more significant dilution of 
the optically thick material around the accretor and evolution 
of a fast wind from the hot star nearby the orbital plane. The 
latter is a direct consequence of the former event, because the 
flux of the hot radiation and particles from the accretor was 
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less blocked in directions of the orbital plane. These events are 
supported by the following results: 

(i) A significant increase in the temperature of the stellar 
component of radiation, Th, we indicated through the model 
SED (Sect. 3.3.3) and by stronger fluxes of highly ionized ele- 
ments (e.g. Fig. 5). Consequently, we have 



(6) 



in contrast to the evolution at/around the maximum (Eq. 5). 

(ii) The optical rebrightening during 2002.7-^2003.3, which 
led to an enhancement of the nebular emission (Sect. 3.6). 

(iii) An additional extension of the Ha wings on our spec- 
trum from 02/02/03. Its temporarily larger terminal velocity 
and flux (Fig. 4) are consistent with the evolution of the hot 
star wind nearby the orbital plane and the enhancement of 
EM (Sect. 3.6). Similar broadening was measured also for 
He 11 ^4686 Hne by Tomov et al. (2005). 

(iv) Appearance of a fast wind from the hot star nearby the 
equator plane was also suggested by the [Fe vii] /i6087 line pro- 
file, which evolved to that arising from the wind-wind interac- 
tion zone (Sect. 3.4.3). 

(v) A sudden increase of the Raman conversion efliciency 
(Sect. 3.4.1, Fig. 8) and the appearance of the pronounced vi- 
olet shifted component in the Raman line profile (Sect. 3.4.2, 
Fig. 9). Both events need an increase of the scattering acts; the 
latter one also negative velocities between the Hi atoms and 
O VI photons. This can be satisfied if a fast wind from the ac- 
cretor and ionizing photons appear near the orbital plane and 
thus interact with the densest portion of the neutral wind from 
the giant around the binary axis that moves against the Ovi 
photons. 

Finally, the model SED from 2003 suggests that we still 
could not see just the bare white dwarf, because of a large ef- 
fective hot steflar radius (Rf ~ 0.1 Rq, Table 2). 

5. Conclusions 

We studied the symbiotic prototype ZAnd by using the far- 
UV, optical low- and high-resolution spectroscopy and UBVR 
broad-band photometry carried out along its major 2000-03 
outburst. We summarize the main results as follows: 

(i) An expansion of the hot object and an enhanced mass 
loss from the system was indicated by the change in the SED 
during the transition from the second plateau stage to the maxi- 
mum (Sect. 3.3.2). A complex mass-outflow from the system 
was recognized in the structure of the [O iii] /15007 profiles 
(ejection of blobs, creation of adjacent emission components; 
Sect. 3.2, Fig. 2) and of the absorption components in C m and 
Pv lines in the FUSE spectra. 

(ii) Evolution of the two-temperature type of the spectrum 
at the initial stages of the outburst having signatures of a mass- 
outflow at moderate and high velocities suggested that the ac- 
tive object consists of an optically thick slowly expanding disk- 
like structured material encompassing the white dwarf at the 
orbital plane and of a fast optically thin wind escaping the star 
at higher latitudes (Sec. 4.1). The optically thick disk-like shell 
persisted probably until 2002 August. Afterwards we indicated 



its dilution and evolution of a fast hot star wind concentrated 
more at the orbital plane (Sect. 4.2). 

(iii) A striking similarity of [Fe vii] /16087 and Raman 
/i6825 profiles at/after the dilution of the disk suggested their 
origin within the interaction zone, where the winds from both 
the stars collide (Sect. 3.4.3). This unusual symbiosis of these 
lines and properties of the Raman scattering process can be 
helping to study streaming along the zone of interaction be- 
tween the two winds. 
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Appendix A: l-ie i flux ratios 

To predict flux ratios of He i emission line s, F(7065)/F(5876) 
andF(6678)/F(5876),Proga et al. ( 1994) solved the statistical 
equilibrium equation for a given level in their (19 level) Hei 
atom including the radiative and three-body recombination and 
eff'ects of collisional ionization and collisional transfer between 
other levels of the model atom. The model requires as input 
parameters the electron temperature, T^, the electron density, 
A^e, the rate of photons capable of ionizing He i atom and the 
optical depth in the He 1/13889 line, T3889, as a measure of the 
population of the metastable 2^ S level. Modeled and observed 
ratios are shown in Fig. A. 1 . The location of the flux ratios from 
the optical maximum in 2000 December suggests that the He -H 
emitting region is very dense, A^e ~ 10'^ cm"^, and radiates 
at Te ^ 20 000 K. Comparing the observed luminosity in the 
/15876 line, AncfF{5dn6), with that predicted by its emissivity, 
£5876 we Can determine the eff'ective radius of a spherical 
nebula as 



/?eff(He 



3F(5^16)d 

Ne £5876 



'2\l/3 



(A.l) 



where the emission coefficient £5376 



^1 X 10 ^^ergcm-^ s ' 
jProga et al. 111994 . 



10'^ cm ^ and T3889 



100 



for 
Then 



F(5876) = 73 X 10-'^ erg cm"^ s^^ (Table A. 1), and c/ = 1 .5 kpc 
yield /?eff(He -h) = 5 Rq. 

During the recent major outburst He i fluxes evolved toward 
lower F(7065)/F(5876), whereas during the smaller 1986- 
eruption the fluxes evolved toward lower F(6678)/F(5876) at 
nearly constant F(7065)/F(5876) (Fig. A.l). In general, this 
reflects different properties of the nebula, which develops dur- 
ing outbursts around the hot active object. Qualitatively, during 
strong outbursts the velocity field of the hot star wind can make 
large differences in the density throughout the nebula. As a re- 
sult, contributions from the outer, low-density part of the neb- 
ula will enhance more the /15876 line relatively to other lines 
and th us make lower the considered flux ratios. According to 
ISmitsi 0996) the relative fluxes F(5876)/F(6678)/F(7065) = 
2.519/0.682/0.637 for A^e = 10^ cm-^ and = 20 000 K (open 
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Fig.A.l. Dependencies of Hei 
andF(6678)/F(5876),for A^e = 



flux ratios, F(7065)/F(5876) 
10'' cm"^ (dotted lines), A^e = 
10"^cm"^ (dashed lines) and Ae - lO'^cm""^ (solid lines), 
for Te = lOOOOK, 15 000K and 20 000 K (light, medium and 
heavy hnes, respectively). Crosses indicate r^ggg - 0.01, 0.1, 
1,10, 100, 1000 from the bottom to the top on the solid heavy 
line (Table 5 of Prqga^eLaL 1994). Observed fluxes are from 
Table A.l. Open circle corresponds to a low-density nebula 
(Ae = 10^cm-\ Te = 20 000 K). 



Table A.l. Observed fluxes of He i emission lines 



Date 


JD24... 


F(5876) 


F(6678) 
[ erg cm"- s" 


F(7065) 

'] 


A/Q* 


10/12/00 


51889.4 


76.7 


60.4 


36.9 


A 


11/12/00 


51890.4 


63.7 


56.4 


37.3 


A 


12/12/00 


51891.3 


77.4 


57.2 


33.1 


A 


30/09/82 


45213* 


19 


19 


18 


Q 


01/06/85 


46218* 


28 


19 


25 


A 


02/12/85 


46402* 


58 


33 


56 


A 


20/10/86 


46724* 


39 


26 


38 


A 


01/12/88 


47497* 


15 


15 


16 


Q 



' A/Q: Active/Quie scent phase 
* from lProsa et^ fl994..^ 



circle in Fig. A.l). Therefore the real size of the He H- nebula 
will probably exceed the effective radius derived for a constant 
density distribution. In addition, its shape is not spherical due to 
the presence of the optically thick disk-like matter at the orbital 
plane, which prevents ionization in these directions (Sect. 4). 



Appendix B: Radius of the Ovi zone 

Radius of the Ovi zone, /?(Ovi), can be obtained from the 
equation of the equilibrium between the flux of photons, 
iph(0 v), capable of ionizing O v ions and the number of re- 





Vo = 27.54x lo" s — 












- 
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3 

Temperature [10 K] 

Fig. B.l. Number of photons capable of ionizing O v ions as a 
function of temperature scaled to the luminosity of 1 000 Lg- 



combinations by O vi ions. For a spherically symmetric nebular 
medium around the central ionizing source we can write 



i'ph(Ov) = y7;(Ovi)3A(Ovi)A+AeaB(Ovi), 



(B.l) 



where A(Ovi) is the abundance of Ovi ions, q'b(Ovi) is the 
recombination coefficient of a free electron with the O vi ion 
for the case B and N+ and Ae are concentrations of pro- 
tons and electrons. In determining Lph(Ov) values we inte- 
grated the Plank curve from the ionization limit of 1 13.896 eV 
(v() = 27.54 X 10'^ s"') to oo. Figure B.l shows the Lph(Ov) 
quantity as a function of the temperature for the luminos- 
ity of the ionizing source of lOOOL©. Further we assumed 
that all oxygen atoms are ionized to Ovi within the zone, 
i.e. A(Ovi) - A(Oi), for which we adopted the solar pho- 
tospheric abundance of 4.6 x 10 -4 (»Asplund et al. "2004') and 
gR(q VI) = 9.5 X 10 '^ cm+^ s'^ for ^ 30 00 K (Gurzadvani 
Il997h . According to ISokoloski et al. I ( l2005l) . the temperature 
of the ionizing source T^'- = 95 000 K and 1 60 000 K during the 
optical maximum and the 2003 post-outburst stage, which im- 
plies Lph(0 V) = 6.4 X 10"^^ and 2.6x10"*^ photons s"' for lumi- 
nosities of 7000 and 4 400 Lq' respectively (Fig. B.l, Table 2). 
Finally, for the particle concentration of the order of lO'^ cm 
which can be expected in a close vicinity of the ionizing source 
due to the mass loss via the stellar wind at rates of a few 
xlO"^ Mq yr"' (Skopal etal..2006.) . one can derive 



R{0 VI) w a few x Rq 



(B.2) 



during the optical maximum. In a more rigorous approach we 
integrate densities of the hot star wind from a radius r\ above its 
origin, /?„, to the end of the O vi zone (i.e. the R(Q vi) radius). 
Then Eq. (B.l) takes the form 



Lpj,(Ov) = 47raB(Ovi)A(Ovi) 



Ji-i 



.«(Ovi) 



N(rf r^dr. 



(B.3) 



where N(r) = M/Anr^fimnvir) is determined via the continu- 
ous mass loss rate M from the hot star wi th the velocity field 
v(r) = Vco(l - R.^/ry^ dCastor et al. Ill975l) . This allows to ex- 
press the last equation as 



iph(Ov) = 



aB(Ovi)A(Ovi) 



\vool Rv, 



X/, 



(B.4) 
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Fig.B.2. Radius of the Ovi zone as a function of Lph(Ov) 
photons for particle distribution in a form of the stellar wind 
from the ionizing source (Eq. B.4). Horizontal bars indicate 
quantities from the optical maximum (2000 December) and the 
2003 post-outburst stage. 



where the integral 
1 



(2/3-1) 



\ n I \ R(Ovi)l 



1-2/3 



(B.5) 



Solution of Eq. (B.4) is shown in Fig. B.2 for three different 
M. Parameters (3 = 1.7, Voo = 2000 kms"', = \ JRq and 
ri ~ 1 .4 Rvj were estimated from the mode ling the broad Ha 
wings during the maximum jSkopal 1 12006). Also in this case 
the R{0 VI) radius can be small, < 3 Rq, for mass-loss rates 
between 1 xlO"^ and 5x10"' Moyr"' and the observed flux 
of ionizing photons Lph(Ov) = 6.4 x 10"*^ s '. However, for 
Lph(0 v) ^ 10'*"^ - lO'*^ s ' the O vi zone can expand to infinity. 
This is in a qualitative agreement with the observed significant 
and fast increase in the O vi /1 1032, 1038 fluxes from the optical 
maximum (2000 December) to the next FUSE observation in 
2001 July (Table 3), which thus can be understood as a result 
of the increase in to 120 000-150 000 K at approximately 
constant luminosity dSokoloski et al. i2005ii) . 
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